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The induced diffusion of tracers in a bacterial suspension is studied theoretically and 
experimentally at low bacterial concentrations. Considering the swimmer-tracer hydro- 
dynamic interactions at low-Reynolds number and using a kinetic theory approach, it is 
shown that the induced diffusion coefficient is proportional to the swimmer concentra- 
tion, their mean velocity and a coefficient /3, as observed experimentally. The coefficient 
(3 scales as the tracer-swimmer cross section times the mean square displacement pro- 
duced by single scatterings. The displacements depend on the swimmer propulsion forces. 
Considering simple swimmer models (acting on the fluid as two monopoles or as a force 
dipole) it is shown that /3 increases for decreasing swimming efficiencies. Close to solid 
surfaces the swimming efficiency degrades and, consequently, the induced diffusion in- 
crease. Experiments on W wild-type Escherichia coli in a Hele-Shaw cell under buoyant 
conditions are performed to measure the induced diffusion on tracers near surfaces. The 
modification of the suspension pH vary the swimmers' velocity in a wide range allowing 
to extract the /3 coefficient with precision. It is found that the solid surfaces modify 
the induced diffusion: decreasing the confinement height of the cell, /3 increases by a 
factor 4. The theoretical model reproduces this increase although there are quantitative 
differences, probably attributed to the simplicity of the swimmer models. 
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1. Introduction 

An active suspension is a fluid containing autonomous swimmers such a bacteria, algae 
or artificial self-propelled entities. In general, the presence of active swimmers changes 
the mechanical picture usually considered for passive suspension problems. When an or- 
ganism swims, it interacts with the surrounding medium (i. e., the suspending fiuid, other 
swimmers in the fiuid and also the boundaries of the system). Consequently, balances 
of momentum and energy as well as the constitutive transport properties, are deeply 
modified by the momentum sources distributed in the bulk. When a microscopic cell 
swims in a fiuid, usually viscous forces dominates the inertial terms, which is evidenced 
if we compare these two magnitudes by computing the Reynolds number Re = pyj±^ 
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where y is a typical velocity, L a typical length, and p and 77 are the density and the 
dynamic viscosity of the fluid, respectively. Swimmers, such as bacteria, with a typical 
length of 1 /im, and a propulsion velocity V ^ 20 fim/s moving in water (p = 1 g/cm '^ 
and 77 = 10~^g/cm-s), leads to a Reynolds number in the order of 10~^ l Purcelllll977t ). 
Under these conditions, au tonomic propulsion is assured only if the time reversibility is 
broken during the motion (|Purcelllll977l : iGolestanian fc Aidarill2008( ). 

Since the swimming motion has no inertia, the total force and torque exerted by 
the swimmer on the fluid vanish. Therefore, the leading long range hydrodynamic flow 
produced by a single swimmer should be governed at most by a dipole force. As a conse- 
quence, depending on the dipole polarity, many micro-organisms can be classified under 
two categori es: "pushers" and "pullers" , depending if the dipole is tensile or contractile, 
respe ctively ([Hernandez-Ortiz et al ISaintillan fc ShellevllioOTilBaskaran fc Marchettil 
From this point of view, Escherichia coli or Bacillus suhtilis can be classified as 



pushers and algae such as Chlamydomonas reinhardtii as pullers. 

Escherichia coli (E. coli) represents a good example of self-propelled swimmer. This 
widely studied bacterium will represent the active component in the solution studied in 
this work. E. coli is a cylindrica l cell with hemispherical endcaps, with 1 pm diameter and 
2 pm length ( Berd 200061 . 2004). The motion is the consequence of the rotation of several 
helicoidal flagella, where each flagcllum is linked to the cell membrane by a nanoscale 
motor. When all the motors turn in a counterclockwise (CCW) direction, the flagella 
rotates in a bundle and this pushes the cell steadily forward (what is called "run"). When 
the motor rotation switches to clockwise (CW), the cell can change direction (known as 
a "tumble" process). The mean run interval is typically 1 s, whereas the mean tumble 
time is around 0.1 s . The ratio between tumble and nm depends upon the chemotactic 
signals (lBerg||2000a|). Due to the drag anisotropy on the helicoidal flagella, its rotation 
produces a net force t hat pushes the fluid backward, while the head pushes the fluid in 
the opposite direction ( Purcelll 1977 : Lauga fc Powers! 20091 ) . As consequence of the fluid 
incompressibility, the fluid is pulled toward the centre of the E. coli body, creating two 
toroidal vortices around it. 

The induced motion of the fluid produced by swimming microorganisms can be mea- 
sured indirectly by observing the induced motion of tracers dispersed in the fluid. The 
seminal work addressing the effect of active suspension on the enhancement of passive 
tracers diffusivity, was presented by Wu fc Libchaber ( 2000l ). They conducted an inter- 
esting experiment: wild- ty pe E. coli as well as 4 .5 and 10 pm diameter particles were 
trapped in a thin film fsee lWu fc Libchabeij[2000l . fi gurcs 2 and 3). When looking at the 
passive tracers, they found the presence of superdiffusion for t < tc and normal diffusion 
for t > tc, where tc is a characteristic time representing the lifetime of coherent struc- 
tures in the sample. They also studied the influence of the bacterial concentration in the 
diffusion of 10 pm diameter particles, showing the diffusion of these particles increases 
linearly with the bacterial concentration. They suggested that this effect could be re- 
lated to the spontaneous formation of swirls in the bacterial bath. In 2009, Leptos et al. 
studied the enhancement of tracers diffusion into an aqueous medium, for a puller type 
swimmer. They used Chlamydomonas reinhardtii as the active component of the solution 
and 2 pm diameter beads as passive tracers. Their observation was performed far from 
the wall and they show a linear increase of the tracers diffusion with the volume fraction 
(/>. They found that the effective diffusivity Dcff for the passive tracers can be written 
as DeS = Db + CK(j), where Db is the Brownian diffusion value in the bulk (without 
swimmers) and, a can be defined (from dimensional analysis) as a C/^t ~ Ul (where 
U, T and / represent a characteristic advective ve locity, encoun t er tim e, and advective 
length, respectively). In a recent work presented by lWilson et all ( 2011 ). Differential Dy- 
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namic Microscopy (DDM) was applied to characterised the bacterial motion in E. coli 
suspensions. They show that the diffusion of non-motile bacteria is enhanced with the 
active fraction, given by the proportion of active swimmers, also establishing a linear 
relationship between these two quantities. 

In 2011, we presented results obtained for passive particles ne ar a solid surface in a 



bacterial suspension of density n under non-buoyant conditions ()Miho et a/.ll2011l ). We 
were able to distinguish two kind of bacterial motion: bacteria that exert a random motion 
and those that move in ballistic way, called active swimmers. From this separation of 
populations sorted according to their motion, one can define a fraction of active swimmers 
(pA^ and their mean velocity Va- To characterise the activity at the surface we define the 
"active flux" J A = tiaVaj where ua = tk^a is the active concentration. Systematic 
experiments varying the bacteria concentration were done for two tracer sizes and in 
each case the mean diffusion of the passive tracer Dp was obtained. Interestingly all 
data seem to group on a straight line 

Dp = D^ + I3Ja, (1.1) 

where Dp is the Brownian diffusivity close the wall in the absence of bacteria. The 
parameter (3 represents the slope that can be determined experimentally by a linear fit. 
A dimensional analysis of the equation (|l.ip yields that this pre-factor is a length to the 
power 4. 

In the present work wc present a kinetic theory analysis considering the swimmer- 
tracer hydrodynamic interactions, that explains the observed dependence of Dp on Va 
and Ua, and provide an expression for /3. When a tracer is placed in the suspension, the 
swimmer-tracer scattering produce a net displacement that, as a consequence of the low- 
Reynolds number hydrodynamics as will be demonstrated in Sec. [H is independent on 
the swimmer velocity. In kinetic theory, the diffusion coefficient can be computed as the 
collision frequency times the net displacement squared. On the other hand, the collision 
frequency is proportional to the swimmer concentration, their mean speed and the cross 
section. Putting all together, results in an expression like (|1.1|) . The /3 coefficient scales, 
then, as the cross section times the average displacement squared. As a result of the 
analysis, it turns out that an enhancement of the induced diffusion can be obtained near 
solid surfaces, calculations that are performed in Sec. [3l Controlled experiments near 
solid surfaces, in which bacterial batches with different swimming speeds are prepared, 
are described in Sec. |4] showing the predicted enhancement of the diffusion near the 
solid surface and under confinement. The numerical comparison show good qualitative 
agreement. Finally, conclusions are presented in SecO 



2. Induced diffusion in the bulk 

Self-propelled objects agitate the fluid when they swim and, therefore, a passive object 
placed nearby will be transported by the fluid. At low Reynolds number, the transport 
is passive with no inertial delay. The experiments described in the Introduction and the 
new ones described in Sec. [4] show that, at low volume fraction, the induced diffusion is 
proportional to the swimmer concentration. This result implies that it should be possible 
to understand the induced diffusion as the additive effect of the interaction of the tracer 
with individual swimmers. Swimmer-swimmer interactions or the collective motion of the 
swimmers would lead to higher order contributions in the concentration to the induced 
diffusion. The hyd rodynamic int e ractio n of passive tracers with a single swimmer was 
studied in detail bv lDunkel eFall (|2010f ). They showed that in the scattering process, the 
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tracer trajectories are loop-shaped, almost closed, with small total displacements. Here 
we consider the effect of these scattering processes in the induced diffusion of the tracers. 

Two simple swimmer models will be considered in order to compute the induced dif- 
fusion on passive objects. At low Reynolds number, swimmers have no inertia and the 
net force and torque acting on them vanish. As a first approximation, a swimmer can 
be modelled as a force dipole. The second, more refined, model consists of two force 
monopolcs separated by a finite distance, representing the head and tail of the swimmer. 

The Stokes equations allow to obtain the fiuid velocity at a point r produced by a force 
monopole F located at r' 

u.,{v) = ,h,{v' -v)F,, (2.1) 
where i and j indicate the Cartesian coordinates, J is the Osscn tensor 

1 c . X 1 1 / „ r, r 



^.(r)^^^.(r)^— + (2.2) 

5ij is the Kronecker symbol, 7 7 is the fluid viscosity and the Einstein summat ion notation 
is used throughout the text ( Happel fc Brenneilll965 : iKim fc Kariliall200fJ ). Thanks to 



the linearity of the Stokes equations, the velocity field produced by a force distribution 
is simple obtained by summing the different contributions to the fiuid velocity. Then, 
two monopolcs located at yb ± an/2 (that is, centred at and separated by a in the 
direction n) produce a velocity field 

u,(r) ^ [J^J{VB + a a/2 - r) - J^j{vb ~an/2~ r)] F^, (2.3) 

that in the dipolar limit reduces to 

Ui(r) = Ji],k{rB - v)aFjnu, (2.4) 

where Jij,k is the gradient of the Oseen tensor in the direction k 

When a passive tracer is placed in the vicinity of a swimmer, it will move with a 
veloc ity equals to the fluid velocity at its centre, corrected by the Faxen contribu- 



tion ([Happel fc Brenneij Il965t iKim fc Karillal 120051) . This Faxen correction is propor 



tional to the tracer r adius squared a nd the Laplacian of the fluid velocity. In the experi- 



ments reported in iMiho et al\ (|201ll ) the induced diffusivity does not present dependence 
on the tracer size for particle diameters of Ifxm. and 2/im, showing that in this regime 
the particle displacements are dominated by the fluid velocity at their center and the 
Faxen correction is subdominant. Accordingly, in what follows we will consider point-like 
tracers and the Faxen correction will be neglected. If precise calculations were needed, 
this correction should be added when the tracer and the swimmer are close, in which 
case the velocity field changes rapidly, but if this would be the case, also a more refined 
swimmer model would be necessary. Then, the instantaneous tracer velocity is the fluid 
velocity at its position given either by (|2.3p or (|2.4p . The tracer motion can be made 
non-dimensional by the following choice of dimensions: lengths are scaled by the swim- 
mer size a,r = r/a, and dimensionless time t is defined in terms of the swimmer velocity 
Vb by t = ta/VB- The resulting tracer equation of motion is 

5 = 7-^G,(?B-?), (2.6) 
at 
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7 = SirVBai]/ F 
is a measure of the swimming efficiency and 



G 



monopolcs / 



g: 



dipolc 



(x) [J,j-(x + n/2) - Jy(x - n/2) 
(x) = J^j,k{:x.)njnk. 



(2.7) 

(2.8) 
(2.9) 



We have assumed, as it is usually the case, that the direction of the force is parallel to 
n, with F = Fn. Therefore, F and 7 arc positive (negative) for pushers (pullers). 

It is worth noticing that the swimming efficiency 7 depends on the geometry of the 
swimmer but it docs not depend directly on the fluid viscosity or the activity of the 
swimmer (e.g. the flagcUa rotation speed for E. coli). Indeed, if the swimmer increases 
its activity, by the linearity of the Stokes equations, it will swim more rapidly keeping 7 
unchanged. Also, if the activity is fixed, and the swimmer is placed in a more viscous fiuid, 
its swimming velocity will decrease and 7 is not modified. In the case of E. coli in water, 
using the typical velocity o f Vb = 20 /im/s (|Bergll2004) . the measured dipo le distance a = 
2 /xm ( Drescher et al. 2011 ). and a fitted value F = 1.13 pN, calculated bv lDunstan et al 
(|2012l ) using a two-sphere model, it gives an efficiency of je.coH = 0.89. The swimming 
efficiency changes, however, close to fixed solid surfaces where, for a fixed force, the 
swimmer moves slower than in the bulk. This issue will be exploited in Sect.[3]as it will be 
shown to increase the induced diffusion. There is, finally, the possibility that dynamically 
the swimmer changes its efficiency when moving in a more viscous fluid. Indeed, it has 
been reported that increasing the viscosity, some bacteria swim faster, effect that is 
explaine d by a modification of the fla.gella geometry due t o a balance of elastic and vis cous 



stresses ( Shneider fc Doetsch 1974 : Atsumi et al. 1996 : Chattopadhvav et al. 2006[) . To 



simplify the analysis and because we have no evidence that there are conformational 
changes in our experiments, we will assume that the geometry remains fixed. 

To integrate Eq. ()2.6p the swimmer motion must be specified. As the dipolar interac- 
tion decays at large distances, the swimmer motion can be assumed to be in a straight 
line, with constant speed, during all the scattering process. That is, we neglect the circu- 
lar motion bacteria perform near surfaces and, consistently, the torque dipole created by 
the rotation of the head and flagella are neglected as well. Finally, we neglect tumbling 
during the scattering process. As usual when describing scattering processes, the swim- 
mer trajectory is parametrized by the direction of motion n, the dimensionless impact 
parameter b = b/a and the azimuthal angle 6, that define completely the vector b that 
goes from the initial position of the tracer to the point of minimal approach of the swim- 
mer trajectory (see Fig. [1}. Then, the full swimmer trajectory in dimensionless units is 
parametrized by = tn + h, with —00 < t < 00, and the origin has been placed at the 
initial position of the tracer. Then, the tracer equation of motion reduces to 

dt 

Figure [5] shows typi cal trajectorie s of tr acers for both models. They show the loop- 
shape also described in lDunkel et m (|2010[ ) for a different model. The net dimensionless 
displacement 



=^-^G^{tn + h-S^) 



(2.10) 



A,: 



rj(+oo) - ri(-oo) = 7' 



dtG,{tn + h-r{t)) 



(2.11) 



is small, effect that can be understood when doing a formal expansion of the equation of 
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Figure 1. Parameters describing the scattering geometry between a swimmer (solid circle) and 
a tracer (open circle) in the bulk (Left) and close to a solid surface (Right). The swimmer moves 
in the n direction (parallel to the surface when close to it). The vector b points from the initial 
position of the tracer to the position of minimal approach of the straight trajectory and its 
magnitude is h, the impact parameter. Close to a solid surface, the tracer initial height is zt, the 
angle 9 gives the horizontal displacement of the swimmer trajectory with respect to the initial 
position of the tracer and the height of the swimmer is 2;b. 




Figure 2. Tracer trajectory that starts at the origin produced by the scattering with a swimmer 
that moves along the x axis from left to right, above the tracer. The dimensionless impact 
parameter is fe = 2, the swimmer efficiency is 7 = 1 and the swimmer model is that of two 
monopoles (Left) and one dipole (Right). The insets show a zoom near the origin to put in 
evidence the small displacements. In both cases, the tracer moves left and down. 

motion for small tracer displacements 

/OO /"CXD /.fi 

dtG^itn + b) - 7-2 / dti G^.jihn + h) dh Gjihn + b) + 0{j~^). 
-00 J ~QC J —00 

It can be directly verified that the first integral in the expansion cancels both in the case 
of two monopoles or a dipole. bideed, this integral corresponds to using the velocity field 
at a fixed position and, as a result of the head and tail pushing in opposite directions, 
the effect of the head and the tail cancel and the total induced displacement vanishes. 

Then, the first contribution to net displacement is proportional to 7^^. Figure [3] shows 
the net displacement computed by numerical integration of Eq. (|2.10p for the two models 
and two values of & as a function of 7. The dependence A ~ 7"^ is obtained from 7 > 0.5 
and not only in the asymptotic regime of large efficiencies. At smaller values of 7 or small 
impact parameters, the displacements are large and the tracer approaches the swimmer 
to the core of the singularities, making predictions not to be trusted. In those cases, 
more refined swimmer models should be used. Note that the net displacement does not 
depend on the sign of 7 at least in the dominant term. Therefore, there is no difference 
between pushers and pullers as far as the net displacement of tracers is concerned. We 
finally remind that the net displacement depends only on b and 7 and is independent on 
the swimmer force or velocity. 

Now, considering a tracer in a low concentration suspension of swimmers, the different 
scattering processes can be studied separately and the total tracer displacement results 
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Figure 3. Net dimensionless displacements A times 7^ as a function of the swimmer efficiency 7 
for two impact parameters: h = 1 and 6 = 2. Two swimmer models are considerer: two monopoles 
(solid line) and a dipole (dashed line). 



in the vectorial sum of the sequence of displacements. Let's call A*^") the displacement 
vector produced in the n-th encounter, then the mean square displacement up the n-th 
encounter is 

2 



(Ai?2) 



(2.12) 



Assuming that subsequent swimmers move in uncorrelated directions, the cross averages 
vanish, resulting in 



{rn) 



(2.13) 



where 6^™^ and 7^™^ describe the m-th encounter, and the factor appears when sub- 
stituting A ~ a A. 

To compute the averages, we assume that swimmers are described by a distribution 
function / such that f {n, F, ^) d'^n dF dj is the number of swimmers per unit volume 
that swim along the direction n, with a force F, and an efficiency 7. We assume that, in 
general, the swimmers in the suspension can have different efficiencies and can swim also 
at different velocities. The tracer induced diffusion coefficient equals one sixth the time 
derivative of the mean square displacement, averaged over the incoming flux of swimmers 
y/(ri, i^, 7) and the impact parameters 



D = ~{AR^) = ^ J A\b,^)f{n,F,^)Vbdbded^ndFd^ 
= Y ^^{b,l)f{n,F,j)Vbdbded^ndFdj. 



(2.14) 



If all swimixiers have the same efficiency 70 the previous expression simplifies to 



(2.15) 
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with 

J A^(b,^o)bdbde, (2.16) 

where n = J f{n,F,'y)d'^ndFd'y is the swimmer concentration and (V) is the average 
velocity in the swimmer sample. The coefficient /3 can be written as 

/3 = (A2)a/6, (2.17) 

proportional to the product of the cross section a and the average square displacement 
over this cross section, which has been shown to be typically small. Therefore fi-^^^ should 
not be understood as a single length scale. 

As has been noted above, the net displacement is roughly proportional to and 
therefore the diffusion coefficient is proportional to 7~^. More inefficient swimmers pro- 
duce larger displacements and, therefore, larger diffusion coefficients. The rationale be- 
hind this is that at low Reynolds number, the flow field and the induced displacement are 
proportional to the swimmer force. On the other hand, the net displacement is propor- 
tional to the total interaction time, that is inversely proportional to the swimmer speed. 
Therefore, a slower swimmer, exerting the same force will induce larger displacements 
because the interaction time increases. 

To compute numerically the coefficient (3 a cutoff distance, 6min should be included 
in the impact parameter: physically, swimmers have a transverse radius that impose a 
minimal impact parameter and, in the models (|2.8p and (j2.9p , a minimal distance should 
be included to not approach the singularity core. Then, A(6) is computed by numerically 
integrating the differential equation (j2.10p and the coefficient f3 is then obtained from 
Eq. (|2.16l) . As the displacements arc larger for smaller impact parameter, the final value 
of /? depends singnificantly on femin and a sensible election should be made. 

In the case of E. coli, 7 = 0.89, and the tracers do not reach the singularity cores 
for impact parameters larger than b = 0.75 in the case of the swimmer modelled as two 
monopoles and b = 1.0 in the case of the dipolc model. These values, which are close to the 
bacterial radius, arc chosen as the cutoff impact parameters. The resulting coefficients are 
^ Q Q3Q ^ 0.48 /xm)-* and /jmonopoic^a" = 0.088 (/3m°n°P°i°^ = 1.40 ^m)^, 
where a = 2 fim has been used. 

As it is discussed in Sec. 21 under our experimental methodology is not possible to 
extract a value for (3 in the bulk to compare with the theoretical predictions. On the other 



hand, analysis of the experimental results on Chlamydomonas reinhardtii bv lLeptos et ah 
( 20091 ) allows to estimate ^C.r. bulk _ j^g g _|_ 0.7 /im)^. This value, however, cannot be 



compared with the predictions of the simple swimmer models described here, that do 
not consider the time dependence of the flow produced by this alga and the higher 
multipole moment s needed to accurately describe the near flow ( Drescher et al. 2010t 
Guasto et alll2010() . 



3. Effect of a solid boundary 

Near a solid boundary, swimmers decrease their swimming efficiency due to hydro- 
dynamic hindrance. Lubrication forces and far-field images modify the drag coefficient 
parallel to the wall such that the swimming speed is 

Vb = a{zB)VB^, (3.1) 

where Vbo is the swimmer velocity in the bulk, is the distance of the swimmer to the 
surface and a is the parallel drag correction factor, that vanishes when the swimmer gets 
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in contact with the surface and goes asympto t ically to 1 at large dista nces ( Brenneijll961 



Goldman et aZ.lll967t iHolmqvist et al\\20od iHuang fc Breuerj|2007t ). According to Eq. 
(|2.7p . the swimmer efficiency is modified by the drag with the wall as j{zb) = ctizB)jo, 
where 70 is the efficiency in the bulk. Therefore, close to a solid surface, the diffusion 
coefficient should increase by a factor a"*. 

There are, however, other effects that should be considered close to a surface that 
can either increase or decrease the diffusion coefficient. First, it is known that E.coli 
and other pusher swimmers a re attracted to and trapped by solid surfaces, and swim 
parallel to them for long time dRamia et all 1993 : Erymier et al. 19951 Lauga et al. 200d 
Berke et a/.ll2008l : Li et al. 2008 : Dunstan et al. 2012). Therefore, we consider trajectories 
that are parallel to the surface at a fixed height zb and with planar orientation n. Also, 
the experimental setups in Hele-Shaw geometries allow to measure the in-plane diffusion 
coefficient and, therefore, only the parallel components of A should be considered. 

The second aspect to be taken into account is geometrical: close to a surface, the cross 
section is reduced as swimmers must be above the surface. In the limiting case of a tracer 
in contact with the surface, the cross section is reduced by one half. We also take into 
account that there is a finite transverse swimmer radius ros, such that if zt is the initial 
height of the tracer, the swimmer height is zb ^ zt + bcos9 (see Fig.[T|). The condition 
Zb > tqb limits the integration region in the b-6 space in Eq. (j2.14p . 

Finally, the presence of the planar s urface modifies th e flow profile induced by the force 
distribution. Using the image method, Blake fc Chwand ( 1974 ) obtained the solution for a 
force monopole, taking into account the non-slip boundary condition, therefore modifying 
the bulk solution ()2.ip . First, the flow field is smaller because momentum is dissipated by 
the surface and, second, the flow field is non-isotropic anymore and depends explicitly on 
the distance of the source to the solid surface. The fiow produced by a force F evaluated 
at ri = (a;i,?/i,zi), is 



Ui(x) = 



1 





















d \ hv ■ 



(3.2) 



where r2 = {x2, 2/2, h} is the position of the force point, r = {xi ~X2, ?/i — 2/2, zi ~h} is the 
vector from the singularity to the observation point, and r = {xi ~ X2,yi — y2, zi + h} is 
the vector from the position of the image to the observation point. In the coordinate axes 
we are using, x and y are the planar directions and z is the perpendicular direction to 
the plane. This expression, simplified after using that the force is parallel to the surface, 
allows to obtain the tracer equation for the dipole and two-monopolc models. 

To study the consequence of these additional effects due to the surface, the (3 coefficient 
will be first computed without considering the modification of the swimmer efficiency 
close to the surface. Again, the dominant contribution to the net displacement goes as 
7~^, but the numerical integration of the tracer equation of motion gives two new features. 
First, if the tracer is close to the surface, the non-slip boundary condition included in 
Eq. p.2p implies small velocities and therefore small displacements. But, the presence 
of the surface breaks the symmetry of the flow and, as a result, the trajectories depart 
more from closed loops as the back-and-forth motion is not symmetric. As a result, the 
net displacement can increase at a flnite distance from the wall and then decreases to 
approach the bulk values. Integrating over the allowed impact parameters and azimuthal 
angles the /3 coefficient is obtained, which shows a maximum for a finite tracer height (see 
FigureS]). This unexpected result, obtained even in absence of the enhanced displacement 




Ql , \ , 1 , \ , 1 

5 10 15 20 
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Figure 4. Effect of the tracer height, zt, on the /3 value. Two swimmer models are considered: 
two monopoles and a dipole. The calculations are done considering the enhanced drag coefficient 
(solid line) or without it (dashed line). The numerical parameters are those that model the E. 
coll, that is, a = 2 fim, 7 — 0.89, tob/cl = 0.25, and the cutoff impact parameter is fomin/a ~ 0.75 
for the two-monopole model and 6min/a = 1-0 for the dipole model . 



due to the swimmer inefficiency near walls shows that the perturbation of the loop-like 
trajectories compensates the other two effects near walls that tend to decrease the tracer 
displacement. 

Finally, when the effect of the enhanced drag coefficient is included, the /3 coefficient 
also shows, with larger numerical values, first an increase of /? with height, reaching a 
maximum value for a finite height, to further decrease for larger distances, approaching 
the bulk value. The results for both swimmer models are presented in Fig.|4]as a function 
of the tracer height zt- 

Considering the previously described numerical values to model the E. coli plus the 
bacterial radius of ros = 0.5 /im, the maximal diffusivity is obtained when the tracer is 
located at a height zt ^ 4.7 /im. At this height (3 is increased by a factor 2 compared to the 
bulk value. That is, the induced diffusivity is doubled close to the surface. Both swimmer 
models, although producing qualitatively similar results, predict different values for /3. 
In this respect, the two-monopolc model [described by the singularity distribution in Eq. 
(|2.8p ] is more precise than the dipolar model as it gives a mor e detailed description of the 
near field (actually, as it was shown bv lDrescher et all (|2010h even three monopoles were 



necessary to describe in detail the near field produced by some microorganisms). Also, the 
two-monopole model allowed to use smaller cutoff values, implying larger displacements. 
In consequence, in what follows, we will concentrate only on the two-monopole model. 

3.1. Confinement by two solid walls 

Now we study the effect on induced diffusion in a confined geometry i.e. when the bacteria 
and the passive tracer are placed between two solid walls. This will correspond to a 
situation studied experimentally in the next chapter. The theoretical calculations are 
similar to the previous one. However they are slightly more involved because the image 
system for point forces between two parallel walls is complex, as the images reflect at 
the walls an infinite number of times, increasing the complexity of the image system 
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Figure 5. Effect of the confining height of the system h on the (3 value. The computations are 
done with the two-monopole model. The main figure shows the predicted values for a buoyant 
tracer averaged over the Boltzmann distribution of tracer heights. The inset shows the predicted 
values for fixed initial tracer heights: zt/o. = 1 (solid line), zr/a = 2.5 (dotted line), and 
ztIo, = 5 (dashed line). The numerical parameters are those that model the E. coli, that is, 
a = 2 fj,m, 7 = 0.89, ros/a = 0.25, and the cutoff impact parameter bmin/a = 0.75. To compute 
the Boltzmann averages, the mean tracer height is that of experiments {ksT /m* g) / a — 1.85. 



at each reflection (jBhattacharva fc Btawzdziewica 120021 ) . As a first approximation, we 
make only one reflection at each wall. That is, the image system presented in Eq. (|3.2p is 
repeated for the top wall as well. The reflection on the top wall breaks the back-and-forth 
symmetry even more resulting in an increase of the net displacements compared to the 
case with just one surface. 

The inset of Figure [5] shows the effect of the confinement height h on the /3 coefficient 
for different fixed tracer heights zt- The main figure shows the effect of confinement 
averaged over the Boltzmann distribution of tracer heights with numerical values given 
by the experiments (details are given in Sec.|4]). First, at short confinement heights, the 
induced diffusivity increases, reaching a maximum and later it decreases to reach the 
value of an unconfined system, subject only to the effect of the bottom surface. Both the 
maximum diffusivity and the height at which it is reached depend on the tracer height. 
In summary, it is predicted then that, close but not in contact to the solid surface, the 
induced diffusivity is enhanced under confinement in Hele-Shaw cells. 



4. Experiments on enhanced diffusion 

To test the theoretical predictions made previously, we performed series of experiments 
on passive tracers activated by swimming bacteria. The goal here is to test the validity 
of Eq. (|1.1[) , varying the bacterial density and the mean velo city over a significant range. 
A crucial point in the previous report bv lMino et all ( 201l[ ) was that a distinction has 
to be made between active swimmers performing ballistic trajectories and the non-active 
ones, displaying diffusive motion. Under these conditions, the effect of enhancement in 
the diffusion of passive tracers can be assessed quantitatively. In the bulk, because of the 
3D motion, bacteria rapidly go out of focus, and there is not enough data to make the 
separation between the two populations. Consequently, we have restricted ourselves to 
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measure close to the surfaces where the separation between active and non-active bacteria 
can be made without ambiguity. Also, in the previous experiments, the dependence of 
the swimmer velocity on the diffusion coefficient was only investigated by taking bacteria 
at different stages of their life cycle i.e. the post division fast swimmers (IN cells) and the 
pre-division slow swimmers (2N cells). In this way, a factor 2 between the corresponding 
mean velocities was obtained. However, IN and 2N populations have different shapes 
which from the previous theoretical analysis, can be a questionable issue since a change 
in shape would imply a change in the efficiency parameter 7. If there were a variability 
in 7 then Eq. (jl.ip is no longer valid and the diffusion coefficient is given instead by the 
more involved expression (|2.14p . Finally, the advantage of the present experiment is that 
by varying the swimming velocity, relation (jl.ip can be studied more extensively. 

Another question addressed in the theoretical sections is the influence of confinement, 
i.e. the presence of another surface limiting the height of the chamber. From the exper- 
iments it is possible to extract the effect of conflnement in the Brownian diffusivity as 
well as in the enhanced diffusion due to the presence of bacteria. 

The bacteria used here are Wild- type E. coli (ATCC 11105). They were grown and 
prepared following the protocol described in Appendix 1X1 

Observation on the bacterial suspension was performed using a Zl inverted microscope 
from Zeiss-Observer. Images and videos were captured using a Pixelink PL-A741-E CCD 
digital camera, connected to and controlled by a computer which stores the images that 
will be post-processed. The CCD chip has a maximum resolution of 1024 x 1280 pixels^ 
and can run at 10 frames/s full frame. To gain speed of acquisition, we reduced the 
visualisation held to 600 x 800 with a rate of 20 frames/s. 

For experiments on induced diffusion, the principle is to follow the motion of passive 
particles in a bacterial suspension. A sketch of the chamber is shown in figu re [71(a). A 10 



lA dr op of the suspension containing the W wild- type E. coli bacteria (see lArcher et 



20111 ) and the tracer latex beads of diameter d = 2fim, were placed in a transparent 



chamber on the visualisation stage of the microscope. The chamber is made of two cover- 
slips separated by a typical distance h, ranging between 5 and 110/zm (see Appendix lAl 
for details). 

As discussed previously our goal is to produce bacterial samples at different con- 
centrations and with different mean swimming velocities. In the work presented by 
Minamino et all (|2003l ) the effect of pH in the motility of E. coli and Samonella motil- 



ity was studied. They varied the pH between 5 and 7.8, showing that in the presence 
of potassium acetate, there is a significant variation of the swimming speed reaching a 
maximum at pH 7. Inspired by this result, we used a similar protocol to vary the velocity 
of the bacteria (see Appendix E| . We obtained a mean velocity variation from 4 /im/s at 
low pH up to 20 /im/s at high pH. that is, an increase of the mean velocity b y a factor 5 
(see F ig. |6l). However it is important to notice that the results presented in lMino et al\ 
( 201lh were obtained in non-buoyant condition using a 1:1 mixture of Minimal M otil- 
ity Medium (MM A) and PercoU, a nanoparticles suspension ( Laurent et al ] ll980h . We 
actually noticed that the viscosity of the MMA-Percoll mixture changes strongly with 
pH and forms a gel at low pH, therefore, we could not obtain a density matching be- 
tween the solution and both the beads and the bacteria while keeping the viscosity fixed. 
Consequently, the experiments presented here were done with MMA only (without Per- 
coU) and we tested by direct rheometric measurements that in all cases the fluid stays 
Newtonian. The pH variation does not affect the viscosity of the fluid with a value, at 
25°C, close to the one of pure water (see Appendix rS| . In this context, the latex particles 
(density 1.03 g/cm'^) will sediment and stay preferentially close to the bottom wall. In 
figure [TJb), bacterial concentration profiles are shown in a situation where the chamber 
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Figure 6. Effect of tlie pH in mean velocity Va of tfie active population using MMA solution. 
Error bars represent standard deviations. 



height is h = 110 /iin in buoyant conditions. We observe an increase of concentration 
close to the bottom surface and the corresponding visualisation region is displayed in a 
rectangle on the same figure. Also an increase of the bacteria concentratio n near the top 



wall c a n be noticed , phen o menon that has been report e d pre viously by iFrvmier et 



(|l995l ): iLauga et al\ (|20Q6l) : ISerke et al\ (|2Q08[) : lLi etai\ (|20Q8l ). as a consequence of the 
trapping of bacteria at the surface by hydrodynamic interactions 

Given the specific protocol explained in Appendix A, once the chamber is placed under 
the microscope, the steady distribution of bacteria in the vertical direction is reached very 
rapidly. The experiments start few minutes later in order to let the latex beads sediment. 
Figure shows a snapshot of several latex beads among the bacteria close to the 
bottom surface and we superpose the corresponding tracks for 30 s. In Fig.[71[d) is shown 
an averaged quadratic displacement as a function of the time lag r that the tracers are 
tracked, obtained at two different frame-rates of acquisition. This indicates that after 
0.5 s , a diffusive motion is reached. The diffusion coefficients extracted by linear fit of 
the quadratic displacement versus time lag are obtained by analysing 300 images, taken 
at 1 frame per second. The time lags used for the fit are between 1-10 s. For all the data, 
the bacterial motion is analysed using 20 s videos at 20 images per second. The active 
bacteria con centration is d etermined by characterising the tracks with the method 
described in lMino etUR (|201ll ). using as indicators the mean persistence angle and the 
exploration radius. For each experiment, the mean bacteria velocity Va is computed from 
the mean track length over the total visualisation time. 

4.1. Diffusivity enhancement close to the bottom surface 

In this first series of experiments the chamber height is maintained at its higher value 
of 110 ^m. As described in detail in Appendix A, the change of pH changes the mean 
bacterial velocity Va- Systematic measurements of the diffusion coefficient of tracers were 
made for suspensions prepared at pH 5.2, 5.5, 5.9, 6.3, 6.9, 7.5, 7.9 and 9. For each pH, 
the bacterial concentration was varied by preparing suspensions at different OD (Optical 
density) ranging from OD = up to OD = 1. As discussed before, for each experiment, 
the corresponding values of the particle diffusivity Dp, the mean active bacterial velocity 
Va and the mean active density ua were measured. Figure [8] shows the relation between 
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Figure 7. (a) Sketch showing a a lateral view of the chamber observed on the microscope. For 
these experiments a 40X objective (Aperture number AN 0.65) was used, given a visualisation 
field of 96 x 128 /im^. Observation are performed in the centre of the drop, far away from the 
drop border. (6) Distribution of bacteria as function of height in buoyancy solution (Minimal 
Motility Medium at pH 7) for two bacterial concentration n: 4.5 x 10*bact/ml (solid squares) and 
7 X 10*bact/ml (empty squares), corresponding to optical densities OD 0.7 and 1, respectively. 
The cell height is /i = 110 pm. (c) picture of tracer trajectories, (d) quadratic displacement as a 
function of time lag at two temporal resolutions. 



passive tracer diffusion Dp and the active flux J a = Jt-aVa. Different symbols represent 
different pH, exce pt for the solid circles which correspond to the results presented by 



Miiio et al\ (120111 ). The data collapse on a straight line as it is proposed in Eq. (|l.ip 



where the Brownian diffusivity near solid boundaries is 

Here, ks is the Boltzmann constant, T is the absolute temperature, 77 the viscosity of 
the liquid, d the diameter of the particle and a < 1 is the parallel drag correction factor 
introduced in p.ip . From the data presented in Fig.[51 the fitted value for a is 0.80±0.08 
and f3 is 13.0 ± 0.7 nin^, where the last one gives the influence of bacterial motion in the 
tracers diffusion. 



Now, comparing with the values reported previously in a non-buoyant condition (jMifio et al. 



201 ih the value a found for passive tracer of 2 fj.m was 0.74 ± 0.03, and for /3 was 
5.0 ± 0.4 /zm^. It can be noticed that for small heights /3 in a buoyant condition it is 
signiflcantly higher than the one found for non-buoyant conditions. This result can be 
observed in Fig. |8l where the previous results in non-buoyant conditions (black circles) 
also are plotted. 
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Figure 8. Enhanced difTusivity Dp of passive tracers as a function of J a in buoyancy conditions. 
Squares represent tracer of 2 diameters. Each symbol represents an experiment performed 
with different pH. Black circles represent the res ult obtained und er isodense (non-buoyant) 
condition for passive tracers of 2 /im diameter (see iMifio 



To make a numerical comparison with the theoretical model, we must note that al- 
though the tracers are buoyant, they are not in contact with the surface due to thermal 
Brownian motion. In fact, the tracer heights are distributed according to the Boltzmann 
distribution p ^ exp{—m*gzT/kBT), where m* is the buoyant mass, with an average 
tracer height kBT/m*g = 3.7 /im. Then, the (3 coefficients shown in Fig. |4]must be av- 
eraged over the tracers distribution, resulting in = 2.0 ^m^ when the two-monopole 
model is used. The predicted numerical values arc smaller than the observations. The 
difference could be attributed to the simplified swimmer models, which do not treat 
adequately the near field flows. 

4.2. Variation of active diffusivity with the confinement 

Now we explore the effect of varying the confinement height h between the top and the 
bottom walls. In these experiments, the pH value is maintained at pH = 6.9 and the OD is 
varied from to 0.9. At each value of /i, we extract by fitting linearly the relation between 
Dp and J a, the Brownian diffusivity coefficient in the absence of bacteria {Dp) and the 
slope P representing the activation effect. The results for various confinement heights h 
are presented in Fig. [31 In the Inset, we plot the relative Brownian diffusivity Dp/Ds, 
where Ds is the theoretical thermal diffusion in the bulk. It can be noted that the 
value of Dp/Ds decreases when the separation between the two wall decreases, which 
is in g ood agreement with the experimental results presented by iLobrv fc Ostrowskv 
( 1996h , in which they measured the transverse Brownian motion of passive particles in a 
confined geometry. In the main figure we display the effect of confinement in the enhanced 
diffusion, given by the slope /3 in Eq. (jl.ip . /? increases significantly when the upper wall 
is near to the bottom one, presenting a maximum at height h ^ 8 fim. This maximum is 
approximately four times the value for unconfined system. 

The observed enhancement in the diffusion under confinement is similar to the pre- 
dicted effect in Sec. 13.11 To make a quantitative comparison, the /3 coefficients showed in 
the inset of Fig.[5]must be averaged over the tracer distribution, results that are displayed 
in Fig. [5l There is a remarkable similarity among the two curves. When decreasing the 
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Figure 9. Effect of the confinement on the j3 value. Experiment was performed in MMA at 
pH 7, using 2 /xm diameter passive particles. Inset: Effect of the confinement in the value a for 
2 /im diameter passive tracers. Dp represent the diffusion of particles without bacteria and Db 
the thermal diffusion in the bulk. 



confining height, both show a maximum in /? at a finite height to further decrease when 
the two plates are very close. The position of the maximum for the model, /i ~ 11 fim, is 
comparable with the experimental value. However, the predicted values and the observed 
enhancement are smaller than the experimental values. This quantitative discrepancy can 
be attributed to the simplicity of the swimmer models, the non-uniformities of the bac- 
terial concentration near the surface and the approximation made to calculate the image 
system between two surfaces (see Sec. I3.ip . 



5. Conclusion and Discussion 

Experimentally, in a mixture of swimmers and tracers, it has been observed that at 
low swimmer concentration, the agitation created by their movement produces an en- 
hanced diffusive motion on the tracers, with a diffusion coefficient that grows linearly 
with the swimmer concentration. A kinetic theory explanation, that uses the properties 
of the hydrodynamic interactions at low Reynolds number, is given to this experimental 
observation. The successive displacements produced by uncorrclated swimmers produce a 
diffusive motion, with a diffusion coefficient that is proportional to the swimmer concen- 
tration, their mean velocity and a constant /3 that scales as the scattering cross section 
and the mean displacement squared in individual swimmer-tracer encounters. A dimen- 
sional analysis shows that the displacements grow for decreasing swimmer efficiencies, 
that is measured as the swimmer velocity times its size and the fluid viscosity, divided 
by the force of the swimmer. It results, then, that inefficient swimmers produce large 
tracer diffusion. The calculations arc done for simple swimmer models described as a 
force dipole or two monopoles. More complex models can be used as well, leading to the 
same qualitative conclusions on the scaling of the diffusion coefficient with the bacterial 
concentration and mean velocity, although with different expressions for /3. A separate 
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analysis, though, should be used in the case of swimmers that have time dependent 
strokes as is the case of Chlamydomonas reinhardtii. 

The swimming efficiency is decreased near solid surfaces, where the modification of the 
viscous drag produces smaller swimming speeds. Near surfaces there are, however, other 
effects that lead to either an increase or decrease of the induced diffusion. Numerical 
calculations, considering all these effects, show that the induced diffusion is enhanced 
near solid surfaces and present a maximum at a finite distance form the wall. 

To verify these results, experiments with Escherichia coli were performed in confined 
Hele-Shaw geometries. Changing the solution pH it was possible to modify substantially 
the swimming speed, and together with a variation in the bacterial concentration, it was 
possible to deduce the (3 coefficient from the measurements of the diffusion coefficient. 
Experiments were done under buoyant conditions and therefore surface effects are present. 
To test the predicted enhancement of tracer diffusion near surfaces, the confinement 
height was changed, finding that /3 increases by 4 times when the cell height is only 8 fim. 
There is a good qualitative agreement with the theoretical calculations, that predict, 
however, smaller increases. 

It must be remarked that the swimmer models used to obtain the numerical predictions 
of /3 were very crude (a dipole of forces or two force monopoles) which are expected to 
describe correctly the fluid motion only at large distances. It is then remarkable that the 
numerical prediction give the correct order of magnitudes and the mean features produced 
by the solid surfaces. To obtain quantitative predictions, the same theoretical framework 
can be used but more accurate swimmer models should be considered, including higher 
multipole moments, to compute the tracer displacements. 

We thank A. Honerkamp-Smith for valuable comments. This research is partially sup- 
ported by Fondecyt 1100100, Anillo ACT 127, and ECOS C11E04 grants. J.D. thanks 
the support of a CONICYT grant (Becas Chile). CM. and E.C. acknowledge the finan- 
cial support of the the Pierrc-Gilles de Gennes Foundation and the Sesame Ilc-de-France 
"SMAC". 



Appendix A. Bacteria growth and working suspensions 

In this Appendix, we describe how the samples were prepared. We present the charac- 
terisation of the minimal Motility Medium (MMA) and the bacteria cultivation as well 
as the beads and chamber preparation. 

A.l. Minimal Motility Media (MMA) 

The choice of the ambient medium is always a difficult issue as it may strongly modify 
bacterial behaviour in terms of chemotaxis, ambient ffuid viscosity or response to pH. 
In the liter ature, several pr o tocols were designed to achieve different goals. For example. 



for E. coli, iBerg fc Brownl (jl972l ) defined an specific protocol to study the chemotac- 



tic effects o f amin o-acids in the bacterial response. A different protocol was chosen by 



Lowe et al\ (|1987I ) as they sought to change the viscosity of the suspending medium in 
order to study the rotation of fiagella bundles with Streptococcus. Note that in general, 
the biological conditions impose a very narrow window of parameters and in practice it 
is always difficult to conduct experiments that change one parameter at a time as soon 
as a set of working conditions is established. In the present study, we consider a Mini- 
mal Motility Medium (MMA) as the swimming medium. The use of MMA has several 
advantages. First, it contains a sufficient amount of nutritional elements to preserve the 
bacterial metabolism, but cell division is strongly reduced, thus allowing control over 
density of the bacterial population and limiting the infiuence of chemotaxis. Second, this 
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medium enables adjustment of bacterial velocity by modifying the pH of the solution, 
acting directly on the molecular proton motor that activates the flagella rotation (see 
Minamino et a/.ll2003l) . A third advantage is that for different values of pH, the MMA 



solution remains Newtonian and it keeps physical properties similar to water as far as 
density and viscosity are concerned. 

MMA was prepared as followed: stock solutions (1 M NaH2P04 and IM K2HPO4) 
were mixed to generate phosphate buffers at a given pH (from 5.5 to 8). The stock 
phosphate buffers were diluted to a fi nal concentration of 10 mM in the pr esence of 0.1 
mM K-EDTA, which favours motility (|Adler fc Templeton|[l967tlAdleilll973[) and 20 mM 
sodium lactate as an energy source. In some cases 1% glucose was added to prevent the 
effects of deoxygenation. Potassium acetate at 34 mM was also present in order to obtain 



bacteria with swimming velocities sensitive to pH (jMinamino et a/.l 120031) 



In order to characterise the viscosity of pure MMA solutions at various pH, an Anton 
Paar Physica MCR Rheometer was used. For each solution, the shear rate was varied from 
10 to lOOs"^, and 20 points were measured at 25°C. The procedure was repeated three 
times. The pH variation does not affect the viscosity of the fluid which stays Newtonian 
with a value, at 25°C, close to the one of pure water . 



A. 2. Bacteria cultivation 

Wild- type E. coli (ATCC 11105) were grown overnight at 25° C in rich medium LB 
(Luria broth). The next day, part of this culture was diluted 100 times in LB, allow- 
ing logarithmic growth to restart. Several protocols were performed in order to obtain 
swimming bacteria with maximal velocity. Some authors claim that the maximum ve 



locity is reached in the mid-exponential phase of the growth curve (jStaropoli fc Alon 



ity 

20001) whereas others, and also our observation, found that new born bacteria fr om the 



post exponential phase run faster (|Amsler et aZ.lll993l: |Prii/j fc Matsumuralll997l) . Thus 



once an optical density of 0.6 was reached (measured with spectrometer from Hitachi), 
the bacterial suspension was washed twice by centrifugation (500 rpm during 10 min) in 
MMA to eliminate all LB traces, and finally resuspended in MMA solution (pH 7) for at 
least an hour, allowing division to end, and obtaining a population enriched in IN bac- 
teria, also called "baby cells" . In some cases, the "baby cells" were sorted by low speed 
centrifugation. The optical density of the resulting suspension was measured in order to 
prepare, by serial dilutions, samples of bacteria at different concentrations. Finally, when 
bacterial suspensions were introduced in the microchamber. Polyvinylpyrrolidone (PVP) 
at 0.005% was added to the suspension to avoid stickiness. 



A. 3. Beads preparation 

To study the effect of bacterial activity on the diffusivity of passive tracers, latex beads of 
2 micron diameter (Beckman-Coulter, density p=1.027 g/ml) were added to the suspen- 
sions. These spherical particles come in an aqueous medium containing surfactant that 
was eliminated by 3 successive centrifugations in MMA. Finally they were resuspended 
in MMA containing 0.005% PVP. 

A. 4. Chamber preparation 

The chamber is made of two coverslips separated by a typical distance h. The largest 
distance h = 110 /xm is achieved by using two calibrated coverslips N°0 (0.085 to 0.13 
mm thick) as spacers. For smaller confining heights, we use as spacer calibrated beads 
previously dispersed in the suspension. The action of capillary forces of the outer rim of 
the droplet will hold tight the bottom plate and the top coverslip at the right distance 
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from each others. In this way we can achieve confinement heights ranging between 5 and 

110/XTO. 

Finally, to avoid the stickiness, coverslips were coated with ethanolic PVP (20 mg/ml 
of PVP in pure ethanol at 20 mg/ml) and dried. 
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